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Abstract. Using the Plateau de Bure interferometer, we observed the A3mm absorption lines of CN, HCN and 
HNC from some of the diffuse clouds which lie toward our well-studied sample of compact extragalactic mm- wave 
continuum sources. The column densities of these species all vary by a factor of about fifty and are prominent in 
only a limited subset of the clouds seen in the most ubiquitous species such as OH,HCO + , C2H and C3H2. We 
searched unsuccessfully for CH3CN and N2H + , which are underabundant compared to dark clouds, by factors 
of at least 10 and 100, respectively. The CN-HCN-HNC column densities vary strongly and non-linearly with 
N(HCO + ), for example, which probably best represents their variation with H2 as well. But their abundances are 
very tightly and linearly coupled to each other, varying in fixed proportion, as is the case for OH and HCO + , 
and (only slightly more loosely) for C2H and C3H2. Having measured one, it is hardly necessary to observe 
the two others in this group. We find (N(HNC)/N(HCN)} = 0.21 ± 0.05, (N(CN)/N(HCN)> = 6.8 ± 1. Such a 
small N(HNC)/N(HCN) ratio is typical of warmer gas in darker, denser environments, and is consistent with the 
notion of molecular formation in warmer media. The 6.8:1:0.21 ratio in diffuse gas is very different from TMC- 
1 (6.8:4.5:4.5) where HCN and HNC are relatively much more abundant. It seems likely that the sequence of 
features with increasing column densities of the CN-HCN family or C n H m -family molecules in diffuse gas actually 
represents a series of gas parcels of increasingly higher H2-fraction, number density, and molecular abundances, 
occurring over a relatively narrow interval of total hydrogen column density. Our experiment seems to have caught 
many molecules in the act of turning on - and turning each other on - in the diffuse interstellar medium. 
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1. Introduction. 



The p resence of interstellar CN was shown by McKellar 
194C|) who identified the R (0) line at A387.4n m (see also 



Swings fc Roscnfcld (1937|) and |Adams (194l|) ). The R(l) 
line being undetectable, and taking as an upper limit on 
its equivalent width a value one-third that of R(0), he 
calculated the "effective" or "rotational" temperature of 
interst 



molecules, are too light for their rotational structure to 
have been used in this way. 

Sixty years after the discovery of the CN lines we are 
still finding out what their presence in diffuse clouds actu- 
ally means. Although it has not generally been regarded 
as much of a challenge to come up with rather conven- 
tional chemical schemes which account for the observed 
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abundances of CN ( Federman et al., 1994 ; Van Dishoeck 
fc Black, 1986), these models beg larger questions which 



of such a temperature in a region with so low a den- 
sity of both matter and radiation has any meaning." Of 
course it does have meaning, and CN has long been used 
to derive the temperature of the cosmic microwave back- 
ground at a wavelength of 2.7mm, the contribution to the 
rotational excitation from collisions often being ncgligi- 



are the subject of this series of papers. Earlier we showed 
that the abundances of the closely related species HCN 
and HNC are much larger than expected in diffuse clouds 
seen toward compact extragalactic mm-wave continuum 



sources at galactic latitudes of 2-15° (Lucas & Liszt, 1993 



bly small ( 


Thaddeus, 1972 


; Crane et al., 198E; Black & 


van Dishoeck, 1991 


). Both CH and CH+, the other early 



1994). Given that the proximate source of CN in conven- 
tional models is the breakup of HCN and/or HNC and 
their immediate antecedents, it follows that our under- 
standing of CN is at best accidental (see also jCrawford & 



Williams (1997)) 
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In the first paper in this new series (Lucas fc Liszt 
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2000 ), we discussed the comparative chemistry of species 



H. Liszt and R. Lucas: Comparative Chemistry of Diffuse Clouds II: CN, HCN, HNC, CH 3 CN & N 2 H H 



bd415+3^9 




--1 



Table 1. Background sources and profile rms 
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Fig. 1. CN, HCN and HNC absorption spectra toward 
B0415+379 (3C111) showing the basic hyperfine struc- 
ture in a fit involving three kinematic components (one of 
which is rather weak) for each species. 



in the C„H m -family of hydrocarbons whose most abun- 
dant known members are CH, C2H and C3H2. Here, we 
present results for the cyanogen- or CN-bearing molecules 
CN, HCN and HNC, based on a modest survey of 
A3mm absorption from clouds occulting a subset of our 
wider sample of compact extragalactic continuum sources 
( Lucas fc Liszt, 1996| ). We find a remarkably tight lin- 
ear correlation among all three, such that it is hardly 
necessary to observe the other two once the abundance 
of any one has been found: a similarly tight linear rela- 
tion exists between OH and HCO + (Liszt fc Lucas, 1994; 



Lucas k, Liszt, 1996 ; Liszt fc Lucas, 2000|) and a some- 
what looser but still quite good one exists between C2H 
and C3H2. By contrast, intercomparisons among the dif- 
ferent families (HCN or C 2 H vs. HCO+; HCN vs. C 2 H) 
reveal a high degree of correlation, but sharply non-linear 
behaviour. In general, all species except OH increase in 
abundance more rapidly than linearly with increasing 
N(HCO+) at N(HCO+) ^ 10 12 cm~ 2 . Analogous non- 
linear behaviour in the optical regime is seen in plotting 
N(CN) vs. N(C2) (power-law slope around 1.4, see below) 
or either of these against N(H2), but tight linear relation- 
ships between species are quite lacking in optical absorp- 
tion studies. 

The manner of taking and analyzing the present obser- 
vations is discussed in Section 2. The main results of this 
work are discussed in Section 3, and a discussion of what 
is known of the chemistry of simple CN-bearing molecules 
is presented in Section 4. 
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0.040 
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Fig. 2. HCO+, CN, HCN and HNC absorption spectra 
toward B0355+508. The HNC spectrum, in which the hy- 
perfine structure only widens and does not create (new) 
components, most readily shows which of the HCO + fea- 
tures are present in the CN family. 

2. Observations. 

The mm-wave absorption data discussed here were taken 
at the Plateau de Bure Interferomer during the years 1993 
- 1997. Table 1 shows the list of background sources ob- 
served, along with their galactic coordinates and the rms 
noise levels achieved in the line/continuum ratio, which is 
also the rms error in optical depth in the optically thin 
limit. The mode of observation was identical to that de- 
scribed previously by Lucas fc Liszt (2000] ) and in our 
other prior work. All data used for the present survey 
were taken with a channel separation of 78 kHz and a 
channel width (FWHM) of 140 kHz, corresponding to ve- 
locity intervals which are given just below for the various 
individual species. 

We have never succeeded in detecting emission from 
HCN, the most optically thick of the species discussed 
here, from the clouds occulting any of our extragalactic 
continuum sources. Thus we derive all column densities in 
the limit of no collisional excitation, i.e. T cxc = T cm b. Very 
weak HCN emission (0.02 K) was seen some 30' South of 
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C Oph at the local HCO + emission peak (0.1 K; Liszt 
(1997|)} but CN emission is absent at very low levels to- 
ward the star itself (Crane ct al., 1989; Kopp ct al., 1996| ) 



2.1. CN 

We observed the (1/2-1/2), (3/2-5/2) and (1/2-3/2) hy- 
perfine components of the N=0-1 absorption line near 
113.490GHz, which have LTE line strengths in the pro- 
portion 988:3333:1235. As shown in Figure 1, the lines 
are split by 3-11 MHz and so are separately resolved. 
Assuming equilibrium with the cosmic microwave back- 
ground, the CN column density is related to the inte- 
grated optical depth of the N=0-1 absorption as N(CN) 
= 2.27 x 10 13 cm" 2 / r -idv. 

Results of Gaussian fitting are given in Table Al (see 
Section 2.6). The channel separation and FWHM of the 
data are 0.206 krns" 1 and 0.370 kms" 1 , respectively. 

2.2. HCN 

We observed the F=l-1, F=2-l and F=l-0 hyperfine com- 
ponents of the N=0-1 transition near 88.63GHz, with LTE 
line strengths in the ratio 1:5:3. As shown in Figure 1, 
the lines are split by 1-3 MHz and are separately resolved, 
though often confused due to blending with other features. 
Assuming equilibrium with the cosmic microwave back- 
ground, the HCN column density is related to the optical 
depth of the N=0-1 absorption integrated over all compo- 
nents as N(HCN) = 1.89 x 10 12 cm" 2 / r -id«. Results of 
Gaussian fitting are given in Table A2 (see Section 2.6). 
The channel separation and FWHM of the data are 0.264 
kms -1 and 0.474 kms^ 1 , respectively. 

2.3. HNC 

We observed the F=l-1, F=2-l and F=l-0 hyperfine com- 
ponents of the N=0-1 transition near 90.66GHz, with LTE 
linestrengths in the ratio 1:5:3. As shown in Figure 1, the 
hyperfine splittings of the J=0-1 absorption of HNC are 
smaller than typical linewidths, and so are observed as 
a broadening of the observed components rather than as 
separate kinematic features. Assuming equilibrium with 
the cosmic microwave background, the HNC column den- 
sity is related to the integrated optical depth of the J=0- 
1 absorption as N(HNC) = 1.76 x 10 12 cm" 2 / r _idu. 
Results of Gaussian fitting are given in Table A3 (see 
Section 2.6). The channel separation and FWHM of the 
data are 0.256 krns" 1 and 0.463 krns -1 , respectively. 

2.4. CH 3 CN 

We observed a band centered at 91.9871 GHz, i.e. at the 
frequency of the strongest hyperfine component of the 
J=5-4 transition in the K=0 ladder, toward B0415+379, 
B0528+134, B1730-130, and B2251+158. This band sam- 
pled as well the same J-transition in the K=l ladder, 



whose intensity should be smaller by a factor 0.076 in equi- 
librium with the cosmic microwave background; yet higher 
K-ladder transitions falling in the band should have truly 
negligible population. Although the spectrum is mildly 
suggestive of a detection toward 3C111, it is not statisti- 
cally significant above the 2a level. The spectrum toward 
3C111 had an rms noise in the line/continuum ratio of 
0.004, leading to a 2-sigma upper limit on the integrated 
optical depth of / rdv < 0.014 krns" 1 , or N(CH 3 CN) 
< 2.4 x 10 11 cm" 2 , N(CH 3 CN)/N(HCO+) < 0.021. These 
are the results quoted in Table 3 and discussed below. 
CH3CN is underabundant compared to the composition 
of dark clouds like TMC-1 and L134N. 



2.5. N 2 H+ 

We observed a band around 93.1738 GHz, the location 
of the J=l-Q transitions of N 2 H+, toward B0355+508, 
B0415+379, B0528+134, B1730-130, and B2251+158, but 
detected no absorption. The best limit is for the gas oc- 
culting B0415+379, and this is the datum used in Table 3 
and discussed below. The spectrum toward 3C111 had an 
rms noise in the line/continuum ratio of 0.0025, leading 
to a 2-sigma upper limit on the integrated optical depth 
of J rdv < 0.016 krns" 1 , or N(N 2 H+) < 2.6 x 10 10 cm" 2 , 
N(N 2 H+)/N(HCO+) < 0.0024. By comparison, the other 
lines of sight did not furnish useful data except to confirm 
that we had not somehow missed large quantities of N2H" 1 " 
anywhere. N2H + is much more strongly underabundant in 
diffuse gas than CH 3 CN. 



2.6. Component fitting 

In Tables Al-3 of the Appendix we give the re- 
sults of Gaussian decomposition, done simultaneously to 
whichever subset of the hyperfine structure was actually 
observed (all for HNC and HCN, and part for CN, see 
Figure 1). The optical depth quoted at line center per- 
tains to the strongest hyperfine component but was de- 
rived from a simultaneous fit to multiple kinematic and 
hyperfine components, the latter assumed to appear in the 
LTE ratio, as discussed in Section 2.4. The integrated op- 
tical depth tabulated is the sum over all hyperfine compo- 
nents for a given kinematic feature, either derived directly 
(for HCN and HNC) or (for CN) by a simple scaling to 
account for that (smaller) fraction of the spectrum which 
was not actually observed. 

We always fit the observed line profile (line/continuum 
ratio) using saturated gaussians of the form exp(— To x 
exp(— (v — v ) 2 /2a v 2 )), simultaneously considering all the 
accessible hyperfine structure. Thus the optical depths and 
linewidths which are presented in the Tables are consistent 
with each other and with the degree of saturation which 
occurs in the data. The fitted linewidths are not artifically 
increased by line saturation. 

As seen in Figure 2, the intrinsic kinematic structure 
along a given line of sight can introduce substantial con- 
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Fig. 3. HCN, CN, HNC, HCO+, C 2 H, C 2 and H 2 col- 
umn densities compared at radio and optical wavelengths. 
Top: N(CN) (open) and N(HNC) (filled) vs. N(HCN). 
Middle: N(HCN) vs. N(HCO+) (open) and N(C 2 H) 
(filled). Bottom, N(CN) vs. N(C 2 ) and N(H 2 ) at optical 
and uu-wavelengths, from the 1995 version of the compi- 
lation originally in Federman et ah (1994). 



fusion as to which clouds contain which molecule. In each 
case, we used our HCO + profile as a template and did 
a simultaneous fit, species-by-species, to all of the hyper- 
fine and kinematic components which were observed. The 
FWHM linewidths in Tables A1-A3 of the Appendix all 
include a contribution from the intrinsic broadening in 
the spectrometer which, as noted above, has channel half- 
power widths somewhat wider than the channel separa- 
tion itself. Where the linewidths are plotted (Figure 5), 
this contribution has been subtracted in quadrature from 
the tabulated value. 

3. Systematics 

3.1. Abundance and linewidth comparisons within the 
family of small CN-bearing molecules 

Figure 3 at top shows that N(CN), N(HCN), and 
N(HNC) have the kind of tight linear relationship pre- 
viously demonstrated for OH and HCO + ( [Liszt fc Lucas 



1996| ; Lucas & Liszt, 1996). The slopes of both regres 



sion lines in Figure 3 are indistinguishable from unity 
(1.03 ± 0.08 in either case) with mean ratios (weighted 
by the variance) of (N(HNC)/N(HCN)) = 0.21 ±0.05 and 
(N(CN)/N(HCN)) = 6.8 ± 1.0. 

Figure 5 shows a comparison of the FWHM linewidths 
derived from profile decomposition (see the Tables in 
the Appendix), corrected for finite resolution. The HCN, 
HNC, and HCO + linewidth intercomparisons show no de- 
viation from the condition of equal linewidths, but CN 
seems to have a slight tendency to be narrower than HCN 
and a more marked tendency to be narrower than HCO + . 
CN is typically observed to be narrower than CH in op- 
tical spectra, where the difference is usually interpreted 
as suggesting a spatial s egregration of C N into the denser 
core regions of a cloud ( Crawford, 1995 ). 

In this work, CN was observed with a somewhat better 
spectral resolution, owing to its higher rest frequency and 
the fixed nature of the channel widths in the autocorrela- 
tor. Although it seems unlikely that resolution dominates 
the linewidth comparison in Figure 5 - HCN and HCO + 
are observed at nearly the same frequency and resolution 
- this should certainly be checked. If CN really is narrower 
than HCN, this is extremely interesting because the tight 
correlation of their column densities militates against any 
really large-scale spatial segregation, such as isolating CN 
in a denser core. Instead, there would have to be a mix of 
species and conditions throughout the gas such that line 
of sight average column densities vary in fixed proportion, 
and the CN linewidth would be smaller throughout the 
cloud, or at least down to scales much finer than we can 
discern directly. 

3.2. Abundance comparisons across families of species 
observed at mm-wavelengths 

Figure 3 in the middle shows the variation of HCN 
with respect to HCO + and C 2 H. As we have shown be- 
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Fig. 4. N(H 2 ) (left) and N(CN) plotted against the to- 
tal H-column density N(H) = N(H I) + 2*N(H 2 ) for lines 
of sight with N(CN) measured at optical wavelengths (as 
in the bottom panel of Figure 3). For convenience, data- 
points corresponding to sightlines with optically-measured 
CN column densities have been outlined. Superposed in 
the left-hand panel are theoretical H2-formation curves 
for uniformly-illuminated gas clots of constant total den- 
sity n(H) as given ( cm~ 3 ) at their left-h and edges. N(H) 
and N(H2) are from Savage et al. (1977 ). 



fore, comparisons across chemical families reveal sharply 
non-linear behaviour; although the variations of HCN 
with HCO + and C2H may appear similar, the C2H- 
HCO + relationship is also non-linear, with the same in- 
crease in N(C 2 H) at N(HCO+) w 1 x 10 12 cur 2 (Paper 
I). C2H however has something of a plateau at small 
N(HCO+), where the mean of N(C 2 H)/N(HCO+) is high. 
The N(HCN)/N(C 2 H) ratio varies from (roughly) 0.08 to 
0.3, increasing with column density. 

Taken as unweighted averages over all observed fea- 
tures, we find <N(HCN)/N(HCO+)>= 1.47 ± 0.86, 
<N(HCN)/N(C 2 H)>= 0.155±0.094; in Paper I we quoted 
<N(C 2 H)/N(HCO+)>= 14.5 ± 6.7 for the unweighted 
mean. When the averages are weighted by the fractional 
errors of the column density, favoring higher column den- 
sity sources, we have <N(HCN)/N(HCO+)>= 1.97 ±0.43 
and <N(HCN)/N(C 2 H)>= 0.23 ± 0.57. 

3.3. Insights from optical absorption data 

Figure 3 at bottom shows optically-determined CN, C2, 
and H2 column densities from the updated version of 
the summary compilation of Federman et al. (1994). It 
is lamentable that there are so few lines of sight for 
which both molecular hydrogen and CN have been ob- 
served. In any case, although the relationships in Figure 
3 at the bottom are non-linear (power-law slopes of 1.42 
and 1.31 reading from left to right) they provide a use- 



ful calibration of absolute abundances. For instance, at 
N(CN) = 2.1 x 10 12 cm" 2 , the fitted curve gives N(H 2 ) = 
4.2 x 10 20 cm' 2 . This can be transferred to the fits at the 
top of Figure 3 to infer N(HCN) = 3.1 x 10 11 cm~ 2 and 
then to Figure 3 at the middle for (very) approximately 
N(HCO+) w 7.2 x 10 11 cm' 2 , or N(HCO+)/N(H 2 ) w 
1.7 x 10~ 9 . This is near the single value N(HCO+)/N(H 2 ) 
= 2 x 10~ 9 that we have used for a variety of purposes 



recently (Liszt & Lucas, 2000), most especially to show 
that the observed amounts of CO in diffuse clouds follow 
directly from recombination of this amount of HCO + with 
electrons in standard cloud models. 

As inferred from this train of thought, the HCO + 
abundance is found to be slightly smaller (w 1.5 x 10 -9 ) 
at N(H 2 ) = 1.0 x 10 21 cm~ 2 . But from this comparison 
we gain the new insight that the abrupt increases in the 
abundances of so many species at N(HCO + ) ^ 10 12 cm~ 2 
are actually occurring at N(H2) ~ 0.3 — 0.5 x 10 21 cm -2 . 



3.4. Linear and non-linear functional dependencies 

It is important to understand which parameters and phys- 
ical conditions are varying in Figure 3. For instance, at 
the top there is a factor 50 variation in N(CN), N(HCN) 
and N(HNC), but much of this - perhaps a factor 20-30 
- occurs over a much smaller range in N(HCO + ): a fac- 
tor 5-6. Thus the total column density of gas probably 
changes far less than might otherwise be guessed. At the 
bottom, N(CN) seems to vary by a factor ten or so at 
N(H 2 ) ~ 3 x 10 20 cm~ 2 , and less rapidly, but still faster 
than linearly, at higher N(H2). 

To tie everything together somewhat more tightly, we 
show in Figure 4 the variation of N(CN) and N(H 2) with 
total gas co lumn density N(H) = N(H I) + 2N(H 2 ) ([Savage 



[et al., 19*77 ) for those sightlines in Figure 3 where CN has 
been measured at optical wavelengths. The point is that, 
while N(CN) and N(H 2 ) are strongly coupled at higher 
N(H.2) in Figure 3, neither N(CN) nor N(H2) increases 
with N(H) or extinction along those same directions. Both 
N(CN) and N(H2) vary substantially over a very small 
range of N(H). [] 

The theoretical H2-formation curves superposed in 
Figure 4 at left were constructed by sampling along the 
mean line of sight through a uniformly-illuminate d spher- 
ical gas parc el of constant total density n(H) ( Liszt & 
Lucas, 2000[ ). They suggest that those sightlines with 
larger N(H2) - and so with larger N(CN) - sample denser 
and perhaps darker gas with a higher local fraction of 
molecular hydrogen. This in turn leads to the idea that 
there is some cascade of effects which accounts for the 



In general the effect of extraneous atomic hydrogen is a 
rightward displacement in Figure 4, causing the data to appear 
consistent with Hb-formation at too low a local number density. 
The anomalous point at large N(H) in Figure 4 for the sightline 
to p Oph A is an extreme example. No other sightline in the 
Copernicus sample had such a large N(H), for any value of 
N(H 2 ). 
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strong changes in column density among the trace 
largely driven by the final approach of the molec- 
rogen fraction toward unity, but perhaps having 
cause the change in some other agent, 
be most classical terms, an increase in local num- 
sity would drive the molecular hydrogen fraction 
it relatively constant extinction, increasing the rate 
ion of ancestor molecular ions like H2CN + , which 
;hen recombine faster at higher density and (if the 
1 pressure is even approximately conserved, as is 
the case in our models) somewhat lower tempera- 
helps considerably that both the electron fraction 
abundance of the basic building material (C + ) 
by the C+/C ratio, which remains high 
gas at the moderate density and extinction rel- 
Here. Earlier we showed that even a constant abun- 
relative to H2) of an ancestor ion could give rise to 
1 han-linear increases in the abundance of a product 
via electron-ion recombination in a gas of con- 
nsity, specifically N(CO) varying as N(H2) 2 when 



Table 2. 



Relative abundances 10 8 x N()/N(H 2 ) 



the 

do ninated 
cliffi se 



from gas with constant N(HCO+)/N(H 2 ) ([Liszt fc 
200C| ). 



-fair ily 



elated matter of interest is the extent to which 
molecules we sample may exist in environments which 
1 partly atomic. In Figure 6 we show the varia- 
the H2-fraction with position inside spherical gas 
of density n(H) = 16 and 128 H-nuclei cm" 3 , for 
ues of the column density across the central line 
through the parcel. At the panel to the right in 
this Figure, where the number density is large enough that 
weak HCO + emission might be seen, even low-extinction 
gas is mostly molecular: the lowest column density consid- 
ered here is about that of a standard "Spitzer" HI cloud. 
But at low density, and when small molecular column den- 
sities are sampled, the ambient environment of hydrogen 
gas is likely to be largely atomic. 

4. Chemistry of simple nitrogen-bearing molecules 

4.1. CN, HCN and HNC 

Shown in Table 2 is a comparison across diffuse and dark 
gas of the chemical abundances of the hydrocarbons stud- 
ied in Paper I and species relevant to the present work 
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represents this work, results from Paper I ( for 
species)), unpublished CH data and a scaling 
!5 N(HCO+) 

TMC-1 are from bhishi et al. (1992|) 



"This work 
C n H 
N(OH) 
Results for 
BD-G is C 
For C OpI - 

jCrawford \c Williams (1997| ) and HCN and HNC are normal- 



Oph model G of |Van Dishoeck fc Black (198fi| ) 
C + is from |Cardclli ct al. (1993|), NH is from 



et al. (1993)) and the claim is made (see the discussion 
of prawford fc Williams (19*97 )) that ordinary gas phase 
diffuse cloud chemistry fails by a factor of about 40 to 
reproduce the observed N(NH) and N(CN). If so, the pre- 
dicted amounts of most molecules would be much smaller 
than arc indicated in the last column of Table 2. 

It seems generally accepted that formation of the CN- 
bearing molecules proceeds via the reaction C + + NH — > 
CN++H, followed by a series of hydrogen abstractions be- 
ginning with CN + + H2 — ► HCN + + H, and terminating 
(forming one of the observed species CN, HCN, HNC) 
after one of many possible electron-ion recombinations 
such as H 2 CN+ + e -> HCN+H ([Van Dishoeck fc Black 



(NH, ( JN, HCN, HNC). The tabulation indicates that the 1986| ; |Fedcrman et al., 1994| ; |Crawford fc Williams, 1997| ) 



abundances of HCN and HNC increase very substantially 
in moving from diffuse to dark gas while CN increases in 
abundance only slightly. 

As represented in the last column of Table 2 by 



the work of Van Dishoeck & Black (1986), models of 
moderate-density quiescent gas-phase diffuse cloud chem- 
istry provide approximately the correct amounts of NH 
and C N, but fail to provide the inferred amount of HCN 



However, the particulars are very much open to question 
since the needed gas-phase reactions do not proceed nearly 
fast enough under typical quiescent conditions in diffuse 
clouds. 

Quiescent gas phase chemical models make NH by cos- 
mic ray ionization of N, followed by hydrogen abstraction 
in the reaction of N + and H2 . This and all other gas phase 



NH formation reactions are discounted by Crawford & 



(which ] we derive by scaling N(CN)) by a factor of 10. 
How ever even this degree of agreem ent may be illusory. 
The |Van Dishoeck fc Black (1986|) results were com- 



puted for a very much larger carbon abundance than is 
now accepted (see the entry under £ Oph from Cardclli 



Williams (1997| ) in favor of formation on grains; they also 
argue that the linewidth of NH, which is typical, rules out 
any formation mechanism involving translationally hot 
particles in the gas phase, like ions (N + ) accelerated in 
MHD shocks. 
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HCN 



FWHM [km s ] HCCT 



Fig. 5. Comparison of FWHM linewidths. Left: HNC 
(top) and CN (bottom) vs. HCN. Right: HCN (top) and 
CN (bottom) vs. HCO+. HCN, HNC, and HCO+ have sta- 
tistically indistinguishable widths. CN is noticeably nar- 
rower than HCO + . 




0.01 



1 

Fractional Radius 



Fig. 6. Theoretical H 2 -formation curves for uniformly- 
illuminated spherical gas parcels of constant density n(H) 
= n(HI)+2n(H 2 ) = 16 cm -3 and 128 cm -3 . In each panel 
the fraction of nuclei in H2, f(H2) = 2 n(H2)/n(H), is plot- 
ted as a function of fractional radius for five values of the 
total column density along the line of sight intersecting 
the center of the parcel, N(H) = 0.25, 0.5, 1, 2 and 4 
cnT 



Whether the NH may be formed in gas-phase reactions 
with fast-moving particles is of interest because a similar 
formation route is cited as the likely source of the unex- 
plained high abundances of HCO+ (i.e. C+ + OH -> CO+ 



+ H) and CH+ (via C+ + H 2 ) ( 


Falgarone et al., 1995; 


Hogerheijde et al., 1995; 


Federman et al., 1996; 


Joulain 



et al. , |_1998j) . One of these pathways also seems to be 
needed to explain the observed amounts of CO, as well, 
and they would help greatly in alleviating the apparent 



shortfall of the CN-bearing species (through the reaction 
of C + and NH). If such processes do not contribute to alle- 
viating the apparent shortfall of NH, it seems possible that 
they might not be needed at all to form the CN-bearing 
species. Alternately, if they do function, they would rep- 
resent strong sources and sinks of NH not considered by 
prawford fc Williams (1997| ). 

Whether the lincwidth of NH rules out an enhanced 
gas phase formation rate really awaits a calculation of 
reaction dynamics and subsequent molecular thermaliza- 



tion in the ambient medium (Bucher & Glinski, 1999). 
For molecules which are formed hot - which includes both 
those which are formed from translationally hot reactants 
and those which carry away significant amounts of energy 
from exoergic chemical reactions - only those which react 
quickly will disappear before reaching kinetic equilibrium 
and so exhibit a broader lincwidth. CH + , for instance, re- 
acts rapidly with both H2 and electrons, while NH does 
not. 

4.2. CH 3 CN 

Table 2 shows that CH3CN is underabundant in diffuse 
gas relative to TMC-1, by at least a factor of twenty, or 
somewhat less, 6, if its abundance relative to HCO + is 
considered. Given that easily-observed species like C3H2 
and HCN are underabundant by a similar factor, the limit 
on CH3CN is not that strong and should probably not be 
considered too surprising. Both of the chemical sequences 



leading to CH3CN in dark gas ( Nejad fc Wagenblast 
|1999| ), (C, CH 2 , HCN, H 4 C 2 N+, CH 3 CN) at early times 
or (C 3 H, C 3 N, HC 3 N+, HCN, H 4 C 2 N+, CH 3 CN) later on, 
run through HCN. 

4.3. A/ 2 H+ 

N 2 H+ is one of several molecules (N 2 H+, NH 3 , S0 2 ) 
which are considered to form only very late in the chemi- 
cal evolution of dark clouds (Nejad & Wagenblast, 1999). 
Observationally, NH3 at least is much more widespread 
than N2H + and has been seen toward B0355+508 by Nash 
(1990D and toward 3C111 n our own unpublished work. 



However, unpublished VLA results communicated by A. 
Marscher indicate that reliable NH3 abundances await a 
survey with the upgraded K-band receivers at the VLA, 
and that neither Nash's profiles nor ours (all taken with 
single dishes), are trustworthy. Unfortunately, NH3 col- 
umn densities are also unavailable from Marscher 's VLA 
profiles, because only the (1,1) lines were observed. 

In molecular clouds, N 2 H + forms in appreciable quan- 
tities at late times via H3 + + N2 — » N2H+ + H2, i.e. when 
the ionization fraction has dropped and when CO has been 
depleted from the gas phase, encouraging a higher abun- 
dance of H3" 1 ". It is known for being detectable only in 
darker regions. 

Toward 3C111, N(N 2 H+)/N(HCO+) < 0.002 (2- 
sigma) as compared to 0.25 in TMC-1. This is by far the 
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best limit we have ever been able to set on an undetected 
species. 3C111 is actually seen through a small hole (re- 
gion of lower than average extinction) in an outlying cloud 
in the Taurus cloud complex; a discussion of the viewing 
geometry is given by 



Lucas fc Liszt (199g). 



5. Summary 

We showed that the abundances of CN, HCN and HNC 
vary in fixed proportion 6.8:1:0.21 respectively. This 
is rather different from conditions in TMC-1, where 
their abundances are in the ratio 6.8:4.5:4.5. The frac- 
tional abundances of HCN and HNC increase substan- 
tially between diffuse and dark gas but that of CN 
increases only by factor of a few. Nonetheless, HCN 
and HNC are much more abundant in diffuse gas than 
predicted by models of quiescent low-density gas-phase 
chemistry, with an intrinsic variation over the approx- 
imate range N(HCN)/N(H 2 ) w 1.4 - 5.0 x 10" 9 for 
N(HCO+)/N(H 2 ) = 2.0 x 10" 9 . The latter value is con- 
firmed by comparing our measurements of CN, HCN, and 
HCO + with the lamentably slight body of optical absorp- 
tion data for lines of sight where both CN and H 2 have 
been observed. 

The linewidths of HCN, HNC, and HCO+ are sensibly 
identical, but CN shows a slight tendency to be narrower 
than HCN and a more pronounced tendency to be nar- 
rower than HCO + . If the CN and HCN are strongly cou- 
pled in column density but different kinematically, this 
argues against a spatial segregation of CN into denser re- 
gions of a cloud (core), such as is usually invoked to ex- 
plain the narrower width of CN absorption lines seen in 
the optical. Instead, clouds would have to be a mixture of 
kinematically (and chemically) distinct regions all occur- 
ring many times along any one line of sight. 

In comparing the abundances of species in the CN fam- 
ily with the hydrocarbons discussed in Paper I (principally 
C 2 H and CaH 2 ) or with HCO + , we find strong correlations 
but sharply non-linear variation. Both the CN-species and 
hydrocarbons show a strong increase in column density 
around N(HCO + ) = 10 12 cm~ 2 and all comparisons across 
families show relative abundance changes by factors of 3 or 
more. We note that C 2 H and CaH 2 vary in fixed propor- 
tion with each other, though not as tightly as for HCO + 
and OH, or as seen here. 

For those sightlines where the column densities of CN, 
HI, and H 2 have been measured optically, it appears that 
N(CN) turns on abruptly at N(H 2 ) w 3 x 10 20 cm" 2 and 
increases faster than linearly at larger N(H 2 ), in a regime 
where the abundance of H 2 itself varies widely over a 
very small range in N(H) (or extinction). Models of H 2 - 
formation suggest that the hydrogen is still largely atomic 
in the gas sampled on sightlines with the smallest column 
densities of the trace species like CN. 

Earlier flLiszt fc Lucas, 2000|) , we showed that N(CO) 
oc 7V(H 2 ) 2 (approximately) when CO forms in a gas of 
constant density via recombination of HCO + , and the lat- 
ter is held in fixed proportion (2 x 10~ 9 ) to H 2 . It follows 



that much stronger variations in trace abundances would 
occur if the higher molecular fractions along some sight- 
lines - reflected in higher N(H 2 ) and N(HCO + ) - arise be- 
cause the density is higher locally; the sequence of molec- 
ular abundances revealed in our work must at least partly 
reflect the tendency for higher-density regions to harbor 
higher H 2 -fractions at a given total hydrogen column den- 
sity. Of course this cannot be the whole story because the 
high abundances of the parent ions needed to reproduce 
our observations of their descendants cannot presently be 
explained by any agreed-upon physical scenario. Equally 
important is to understand why the abundances of some 
species relative to H 2 vary so strongly (e.g. C 2 H, C3H 2 , 
CN, HCN, and HCN) while those of CH, OH, and HCO+ 
seem nearly fixed. 

We failed to find CH 3 CN and N 2 H+, meaning that 
they are at least 10 and 100 times less abundant (relative 
to HCO+) in diffuse clouds than in TMC-1. The limit on 
CH3CN may not be very significant, given that CsH 2 and 
HCN are underabundant by about the same amount but 
are readily seen along many lines of sight and that, in dark 
gas, the formation paths to CH3CN pass through HCN at 
both early and late times. The limit for N 2 H + may be 
more interesting but this species is notoriously hard to 
detect even in translucent clouds of much higher density 
than we study. 

We discussed what is known of the chemistry of the 
CN-HCN-HCN species. Although there is a consensus that 
they form through some chemical pathway involving NH, 
most likely interacting with C + , the origin of NH and the 
mode of interaction are entirely uncertain. 

The next papers in this series will discuss the sulfur- 
bearing species CS, SO, S0 2 , H 2 S and HCS + , and proba- 
bly NH 3 and H 2 CO. 
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Table A.l. CN Absorption line decomposition products 



Source 



B0212+735 
B0355+508 



B0415+379 



B0528+134 

B1730-130 

2200+420 



km s 1 
3.612(0.033) 
-17.076(0.024) 
-11.197(0.380) 
-10.295(0.023) 
-8.292(0.073) 
-2.458(0.017) 
-0.872(0.011) 
0.246(0.063) 
9.582(0.030) 
4.980(0.139) 
-0.974(0.038) 
0.019(0.069) 



1.536(0.104) 
0.995(0.073) 
0.115(0.031) 
1.432(0.096) 
0.267(0.041) 
1.431(0.030) 
3.459(0.060) 
0.169(0.027) 
0.220(0.010) 
0.049(0.008) 
0.852(0.056) 
0.367(0.046) 



FWHM 
km s _1 
0.759(0.061) 
0.532(0.045) 
1.621(0.696) 
0.579(0.051) 
0.706(0.159) 
0.951(0.034) 
0.717(0.019) 
0.459(0.142) 
0.788(0.065) 
1.268(0.312) 
0.687(0.080) 
0.635(0.154) 



j rdv 
kms -1 
2.067(0.217) 
0.939(0.106) 
0.331(0.172) 
1.472(0.163) 
0.335(0.092) 
2.415(0.100) 
4.399(0.140) 
0.138(0.048) 
0.308(0.029) 
0.111(0.033) 
1.038(0.139) 
0.413(0.113) 



Table A. 2. HCN Absorption line decomposition products 



T 



Source 



TO 



kms 



FWHM 
kms - 



rdv 
kms - 



B0212+735 
B0212+735 

B0355+508 
B0415+379 



B0528+134 
B1730-130 

B2200+420 

B2251+158 
B0355+508 



-0.092(0.044) 
3.481(0.011) 
-0.104(0.047) 
2.655(0.101) 
3.517(0.019) 
-17.205(0.009) 
-13.770(0.049) 
-10.406(0.006) 
-8.397(0.027) 
-4.632(0.073) 
-3.975(0.274) 
-2.518(0.016) 
-0.926(0.011) 
9.609(0.014) 
4.914(0.070) 
5.983(0.199) 
-1.426(0.018) 
-0.527(0.035) 
-9.545(0.048) 
-17.207(0.007) 
-13.727(0.050) 
-10.407(0.005) 
-8.433(0.025) 
-4.633(0.066) 
-2.134(0.000) 



0.234 
2.584 
0.238 
0.204 
2.750 
1.345 
0.191 
1.828 
0.386 
0.163 
0.107 
2.540 
5.442 
0.411 
0.094 
0.034 
1.111 
0.451 
0.054 
1.390 
0.164 
1.851 
0.367 
0.154 
0.095 



(0.031) 
(0.037) 
(0.031) 
(0.032) 
(0.044) 
(0.036) 
(0.010) 
(0.032) 
(0.011) 
(0.007) 
(0.016) 
(0.045) 
(0.074) 
(0.007) 
(0.006) 
(0.005) 
(0.015) 
(0.014) 
(0.003) 
(0.017) 
(0.007) 
(0.026) 
(0.009) 
(0.006) 
(0.000) 



0.578 
0.890 
0.605 
0.811 
0.805 
0.603 
1.301 
0.547 
1.228 
2.104 
1.867 
0.866 
0.818 
1.045 
0.894 
0.986 
0.785 
0.799 
1.169 
0.583 
1.419 
0.542 
1.335 
2.156 
1.441 



(0.099) 
(0.018) 
(0.104) 
(0.143) 
(0.026) 
(0.015) 
(0.123) 
(0.010) 
(0.051) 
(0.154) 
(0.486) 
(0.032) 
(0.017) 
(0.029) 
(0.128) 
(0.382) 
(0.027) 
(0.053) 
(0.105) 
(0.012) 
(0.125) 
(0.009) 
(0.048) 
(0.139) 
(0.000) 



0.259(0.056) 
4.406(0.107) 
0.276(0.060) 
0.318(0.075) 
4.241(0.154) 
1.555(0.057) 
0.477(0.051) 
1.915(0.049) 
0.907(0.046) 
0.655(0.056) 
0.382(0.115) 
4.216(0.174) 
8.532(0.212) 
0.824(0.027) 
0.161(0.026) 
0.065(0.027) 
1.670(0.063) 
0.690(0.050) 
0.120(0.013) 
1.553(0.036) 
0.445(0.044) 
1.923(0.041) 
0.940(0.040) 
0.637(0.048) 
0.262(0.000) 
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Table A. 3. HNC Absorption line decomposition products 

Source v T FWHM j rdv 

kms -1 kms -1 kms -1 



B0212+735 -0.512(0.035) 0.341(0.069) 0.226(0.048) 0.148(0.044) 

3.766(0.042) 0.455(0.024) 1.159(0.107) 1.009(0.107) 

B0355+508 -17.172(0.004) 0.383(0.006) 0.577(0.014) 0.423(0.012) 

-13.713(0.104) 0.016(0.001) 1.910(0.270) 0.058(0.009) 

-11.229(0.150) 0.020(0.003) 0.735(0.355) 0.029(0.015) 

-10.384(0.005) 0.676(0.009) 0.498(0.014) 0.645(0.020) 

-8.564(0.022) 0.074(0.001) 1.577(0.060) 0.223(0.009) 

-4.241(0.063) 0.024(0.001) 1.715(0.156) 0.080(0.008) 

B0415+379 -2.467(0.004) 0.367(0.002) 1.004(0.011) 0.706(0.009) 

-0.874(0.002) 1.662(0.006) 0.743(0.006) 2.365(0.021) 

0.127(0.021) 0.083(0.005) 0.491(0.064) 0.078(0.011) 

B0528+134 9.382(0.052) 0.105(0.008) 0.964(0.146) 0.193(0.033) 

B1730-130 4.890(0.093) 0.022(0.003) 0.897(0.275) 0.038(0.013) 

6.733(0.305) 0.007(0.003) 0.897(0.918) 0.011(0.013) 

B2200+420 -1.209(0.061) 0.217(0.019) 1.015(0.164) 0.423(0.078) 

B2251+158 -10.092(0.190) 0.018(0.004) 1.305(0.492) 0.046(0.019) 



